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Abstract 

Random 1/2-disk routing in wireless ad-hoc networks is a localized geometric routing scheme in 
which each node chooses the next relay randomly among the nodes within its transmission range and 
in the general direction of the destination. We introduce a notion of convergence for geometric routing 
schemes that not only considers the feasibility of packet delivery through possibly multi-hop relaying, but 
also requires the packet delivery to occur in a finite number of hops. We derive sufficient conditions that 
ensure the asymptotic convergence of the random 1/2-disk routing scheme based on this convergence 
notion, and by modeling the packet distance evolution to the destination as a Markov process, we derive 
bounds on the expected number of hops that each packet traverses to reach its destination. 

Index Terms 

Random Geometric Routing Schemes, Guaranteed Delivery, Convergence, Markov Process. 

I. Introduction 

A wireless ad-hoc network consists of autonomous wireless nodes that collaborate on com- 
municating information in the absence of a fixed infrastructure. Each of the nodes might act 
as a source/destination node or as a relay, where communication occurs between a source- 
destination (S-D) pair through a single-hop transmission if they are close enough, or through 
multi-hop transmissions over intermediate relaying nodes if they are far apart from each other. 
The selection of the relaying nodes along the multi-hop path is governed by the routing scheme. 
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Takagi and Kleinrock (U introduced the first position-based routing scheme, coined as Most 
Forward within Radius (MFR), based on the notion of progress: Given a transmitting node 
S and a destination node Dst, the progress at relay node V is defined as the projection of 
the line segment SV onto the line connecting S and Dst. In MFR, each node forwards the 
packet to the neighbor with the largest progress, or discards the packet if none of its neighbors 
are closer to the destination than itself. Takagi and Kleinrock [1J also discussed a method in 
which one of the nodes with forward (i.e., positive) progress is chosen at random, arguing that 
there is a tradeoff between progress and transmission success. This scheme bears the closest 
resemblance to the routing scheme that we consider in this paper, where each node picks the 
next relay uniformly at random among the nodes in its transmission range over the 1/2-disk 
oriented towards the destination. There are some other variants of MFR in the literature. In 
0, the authors introduced the Nearest Forward Progress (NFP) method that selects the nearest 
neighbor of the transmitter with forward progress; in [15], the Compass Routing (also referred to 
as DIR method) was proposed, where the neighbor closest to the line connecting the sender and 
the destination is chosen. 

All the aforementioned work lacks rigorous analysis of the conditions for guaranteed packet 
delivery and the expected hop-count between the S-D pairs. In particular, most of the existing 
convergence and hop-count results, based on the notion of progress, ignore the vertical distance of 
the next hop to the line connecting the current node to the destination, which has a considerable 
effect on the packet delivery as the packet gets closer to the destination. For other relaying 
schemes, there have been some results on routing convergence analysis. In particular, Xing et al. 
[(4l showed that the delivery can be guaranteed between any S-D pair using the greedy forwarding 
scheme (GFR), if the transmission radius is larger than twice the sensing radius in a fully covered 
homogeneous wireless sensor network. In [,5J, the authors derived the critical transmission radius 
that ensures the network connectivity asymptotically based on the GFR scheme. 

In this paper, we elaborate on the sufficient conditions for asymptotically guaranteed packet 
delivery of the random 1/2-disk routing scheme and derive bounds on the expected number of 
hops that a packet traverses to reach its destination. The remainder of this paper is organized 
as follows. In Section [TT] we introduce the system model and define the notion of convergence 
for geometric routing schemes, based on connectivity and finite hop-count properties. Then, we 
explain the localized properties that ensure the convergence of the random 1/2-disk routing 
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scheme with probability approaching one as the network density goes to infinity. In Section Hm 
we establish conditions on the node transmission range that ensure connectivity for the 1/2-disk 
routing scheme in a dense network. In Section [IV] we derive bounds on the expected number of 
hops between an arbitrary source-destination pair and show that all packets will almost surely 
reach their destinations in finitely many hops if the network is connected. Section |V] concludes 
the paper. 

II. System Model 

Let \A\ be the area of a convex region A C IR 2 over which the network resides. Network 
nodes are distributed according to a homogeneous Poisson point process with density A. Each 
node picks a destination node uniformly at random among all other nodes in the network, and 
operates with a transmission range R. 

For a general geometric routing scheme, when the destination is out of the one-hop transmis- 
sion range of a given source/relay node, the next relay is selected (based on some rule) among 
the nodes contained in the relay selection region (RSR) of the transmitting node, where the 
RSR in general can be any subset of a full disk of radius R centered at the transmitting node. 
We define the rule that governs the selection of the next relay in each node's RSR as the relay 
selection rule (RSL). Note that a node can potentially act as a relay only if it is contained in 
the RSR of the current source/relay. We define the convergence of a routing scheme as follows. 

Definition 1. A routing scheme is convergent if the following two conditions are satisfied: 

1) Connectivity: Between each pair of nodes in the network that are not already in the 
transmission range of one another, there exists a sequence of relays complying with the 
RSL. 

2) Finite hop-count: Given such sequences determined in the connectivity condition above, 
their cardinalitieW are finite. 

Specifically, in this paper we study a special case of the localized geometric routing schemes, 
namely the Random 1/2-disk routing scheme, where for each transmitting node S in the network, 
as illustrated in Fig. CD the next relay V is chosen randomly among the nodes contained in the 



'Measured by the number of hops/relays. 



Fig. 1. The random 1/2-disk routing scheme. 



1/2-disk of radius R centered at node S and oriented towards the packet destination Dst. Hence, 
the name random 1/2-disk routing scheme is coined. We denote the relay selection region and 
the relay selection rule of the random 1/2-disk routing scheme by |RSR and rRSL, respectively. 
Observe that according to our routing scheme, the next chosen relay might be farther away from 
the destination than the current node. This fact differentiates our routing scheme from GFR [@] 
and its variants. 

Note that for the network to be connected based on any geometric routing scheme, the 
connectivity condition in Definition Q] requires the existence of RSL-compliant routing paths 
between all source-destination pairs in the network. In particular, the network is connected 
under the random 1/2-disk routing scheme if the |RSRs of every node S in the network contain 
at least one node with high probability (note that the |RSR of node S can be oriented in any 
direction that Dst may exist). This ensures the existence of a relaying node in every direction 
of a particular transmitting node and ascertains the possibility of packet delivery to a particular 
destination from any direction. In fact, a node has a possible relay/destination in any direction 
if and only if it can serve as relay/destination from any direction. Now, given that the network 
is connected, the packets are guaranteed to be delivered^ to their destinations, according to the 
random 1/2-disk routing scheme, if the cardinality of the rRSL-compliant relay sequence (i.e., 
the number of hops) starting from the source and ending at the destination is finite almost surely. 
In the following, Theorem Q] presents the main result of this paper regarding the convergence 

2 From here on we claim that the following two statements are equivalent: "the packet delivery is guaranteed" and "the routing 
scheme is convergent". 
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of the 1/2-disk routing scheme. For notational convenience, we let N = X\A\ designate the 
expected number of nodes in the network region of area \A\ and d = z p^- denote the normalized 
area of a full disk of radius R relative to the area of the whole region, such that dN is the 
expected number of nodes in such a disk. Also, f(n) = o (g(n)) means that there exist positive 
constants c\ and k such that f(n)/g(ri) < c\ whenever n > k, and f(n) = Q (g(n)) means that 
there exist positive constants c±, c 2 and k such that c\ < f(n)/g(n) < c 2 whenever n > k. 

Theorem 1. If d = o (N~ 2 / 3 ) and dN — > oo as N — > oo, then the random 1/2-disk routing 
scheme is convergent, according to Definition [7] with probability approaching one. 

Proof: Here we only sketch the outline of the proof and present the respective details in 
the following sections. According to Definition [Q we need to show that asymptotically (i.e., 
as N — > oo) a) the network is connected based on the random 1/2-disk routing scheme; and 
b) given a connected network, the random 1/2-disk routing scheme delivers each packet to its 
destination in a finite number of hops. 

In Section HTH regarding condition a) described above, we derive the upper bound <jn on the 
probability that some nodes in the network are not connected to the rest of the network. Then 
we determine conditions that ensure the convergence of a N to zero as N — > oo. This establishes 
the sufficient conditions for the network to be connected with probability approaching one as 
N -> oo. 

In Section [TV] regarding condition b) described above, we show that given a connected 
network, the number of hops that each packet traverses to reach its destination with the proposed 
routing scheme is finite almost surely, by modeling the distance evolution of the packet to its 
destination as a Markov process. This, together with the results in Section [nil concludes the 
proof for the asymptotic convergence of the random 1/2-disk routing scheme. ■ 

III. Sufficient Conditions for Network Connectivity 

As mentioned earlier, the network is connected if there exists a sequence of rRSL-compliant 
relays between any two nodes of the network, which are not already in the transmission range of 
each other. Hence, the problem of proving connectivity amounts to bounding the probability a N 
that there exist some network nodes which are not connected to the rest of the network. A node 
is connected to the rest of the network if it has a possible relay/destination in any direction, 
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which is equivalent to being able to serve as relay/destination from any direction. Thus, a N 
equals the probability that, for some network nodes, there are certain directions at which their 
|RSRs are empty. We can identify two types of network nodes based on their distances to the 
edge of the network: Nodes that are farther than R away from the edge of the network, which 
we call interior nodes, and nodes that are closer than R to the edge of the network, which we 
call edge nodes. 

For interior nodes, it is clear that the node distribution in their |RSRs, looking in any 
direction, is similar. Therefore, the existence probability of an empty |RSR for an interior node 
is independent of its direction. However, due to the proximity of edge nodes to the boundary 
of the network, the existence probability of an empty |RSR for an edge node highly depends 
on its orientation. For example, the |RSRs that fall partly outside the network region are more 
likely to be empty than the ones that are fully contained in the network region. Hence, we derive 
the probabilities of a node being disconnected from the rest of the network separately for the 
interior nodes and the edge nodes, denoted by a' N and a'^, respectively. 

We now slightly generalize the problem as follows. Suppose that the RSR of a node is 
instead a wedge of angle 2%r] and radius R, with < 77 < 1 (hereafter called ??-disk or 
77RSR, interchangeably). Hence, the ^RSR is a special case of 7/-disk with rj = 1/2. Each 
//-disk has an expected number of nodes rjdN. As shown in Section IIII-Cl the probability of 
network disconnectivity increases as r\ decreases. However, we conjecture that the packet reaches 
its destination with fewer hops as r] decreases. Hence, there may exist a tradeoff between the 
probability of network disconnectivity and the expected number of hops between a S-D pair 
parameterized by 77. We leave the study of this tradeoff to a future work and only derive (in 
Section HVl) the expected number of hops between a source and its destination for the case of 
V = 1/2. 

A. Calculation of a' N 

Consider an interior node x, fixed for now. Given i > 1 nodes in the transmission range of 
x, their directions (in reference to x) are independent and uniformly distributed on [0, 2ir]. The 
probability that x is not connected to the rest of the network equals the probability Uiijf) that 
the angle of the widest wedge containing none of these i nodes is at least 2i]7r. It is not difficult 
to give a simple upper bound on Ui{r]): Of the i nodes, without loss of generality (W.L.O.G.), 
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Fig. 2. A realization for which the widest wedge between the nodes is of an angle at least 2r]n. 

we can assume that (at least) one is at one end of an empty wedge with angle of 2?77r, while the 
other i — 1 are distributed independently and uniformly in the remainder of the full transmission 
disk, as shown in Fig. [2] Note that, by assuming uniform (and independent) distribution for 
the remaining i — 1 nodes over the remainder of the transmission disk, we are also taking into 
account the possibility of an effective empty wedge with an angle 2v(tx greater than 2i]n, as 
suggested in Fig. [2] Hence, we obtain Ui(rj) < i(l — rj) l ~ l . Of course, if i = the probability 



One can obtain the a more precise expression for Ui(r]) using results in [0, page 188 as: 



where [^J is the largest integer smaller than x. This expression is based on the inclusion-exclusion 
principle for the probability of the union of events, for which the first summand provides an 
upper bound and the first two summands provide a lower bound. 

Therefore, averaging over i (number of the nodes in the transmission range of x) and over 
the number of interior nodes, we have: 



is U (rj) = 1. 
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(1) 



where (1 — (2 — \fd)\fd)N is the expected number of interior nodes if, for simplicity, we assume 



B. Calculation of a'^ 

So far we have considered the interior nodes which are well away from the boundary of 
the network region. Now, we consider edge nodes which are closer than R to the edge of the 
network and do not have any nodes in the portion of their transmission ranges that falls outside 
of the network region. Therefore, some ?7-disks (of an edge node) may fall partially (up to half) 
outside the region, which increases the chance that they are empty. We refer to this phenomenon 
as edge effect. Assuming the network region is circular, the number of such edge nodes equals 
(2 — y/d)VdN, which is of order 6 (^/dN^j. We need to determine how their contribution to 
the bound on probability of connectivity differs from the previous computations. 

Consider an edge node e, (5'R) -distance away from the network edge, with < 8' < 1. 
Define angle to be perpendicular to the network edge and pointing towards the edge with node 
e as the origin; in other words, as shown in Fig. [3] we take node e as the pole and the ray ev 
(perpendicular to the network edge) as the polar axis of the local (polar) coordinates at node e. 
We observe that, due to the curvature of the network edge, the overlap of node e's transmission 
range with the network region is larger than what it would be if the network's edge were straight 
(i.e., the line passing through the intersection points A and B in Fig. [3]). This area difference 
(the shaded area in Fig. |3]) is no larger than \fdR 2 containing an expected number of nodes on 
the order of (<i 3 / 2 iV), where the maximum area difference is obtained when node e is located 
on the straight network edge approximation (i.e., in the middle of AB in Fig. [3]). Accordingly, 
we make a further simplifying assumption that d = o (A r " 2 ' 3 ) ; this is equivalent to a practical 
assumption that the ratio between the transmission range R and the radius L of the network 



that the network region is a disk of radius L = *J\A\pK = R/y/d. 
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Fig. 3. Edge curvature. 

region goes to zero fast enough such that the expected number of nodes in the shaded area of 
Fig. [3] goes to zero as N — > oo. Then the error in calculating the probability of any event in the 



we proceed as if the network region is straight where it intersects with the node's transmission 
disk for large N. 

We argued in the beginning of this section that, for edge node e, the probability of an //-disk 
being empty, depends highly on its orientation. Let us consider this claim more closely. Let 
(p = cos _1 (<5) G (0, 7r/2), as shown in Fig. HI where 5R is the distance between node e and 
the straight approximation of the network edge as defined before with < 5 < 1. Note that all 
the //-disks are oriented towards the destination of the packets. Hence, for all 7/-disks that are 
oriented toward an angle in the range (—ip, </?), we must have that the destination is within node 
e's transmission range. Therefore, we only need to be concerned with empty 7/-disks oriented 
towards an angle in the range (</?, 2n — </?). The //-disks oriented to an angle in either the range 
(— cp — 7]tt, —(p) or the range ((p, (p + r]7i) are partially outside the network region, as illustrated 
in Fig. |4l and those oriented to any angle in (cp + rjn, 2n — cp — rjn) are fully contained inside the 



following will be a factor of no more than e 



kd s / 2 N 



— > 1, where A; is a finite constants Henceforth, 



3 This rate will apply as long as the region is convex and smooth. 
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network region; this defines the effective 77RSR range for node e. Note that here, all the angles 
are measured relative to the polar axis ev. 

We now compute a'^ for node e. First, suppose that there are no nodes within the transmission 
range of node e; this event occurs with probability no greater than e~ dN ^ 2 , which corresponds 
to the case where node e's transmission disk is half outside the network region. Therefore, the 
probability that some edge nodes having no other nodes within their transmission ranges is 
upper-bounded by 2\fdNe~ dN l 2 . 

Second, suppose that there are i > 1 nodes in the intersection of node e's transmission range 
with the network region. If an empty ?7-disk exists and it is completely contained within the 
network region, then, W.L.O.G., there should be a node on its counter-clockwise edgdj at some 
angle 9 E (if, 2n — tp — 2t]tx). However, for an empty 77-disk that is partially contained within the 
network region there should be, again W.L.O.G., a node at an angle 9 E (—99 — 2?77r, — if — 7771") or 
9 E ((p — r]TT, (f) on the counter-clockwise edge of the 77-disk. Clearly, the existence probability 
of such empty 77-disks (that is partially contained in the region) increases as either 5 or \9\ 
decreases. The area of the intersection between the 77-disk (that is partially contained in the 



4 The straight edge (of the 77-disk) that over its counter-clockwise direction, the entire 77-disk is located, as shown in Fig. [4] 
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region) and the network region is that of a wedge with angle \9\ — if (wedge eAB in Fig. 
IU) plus a triangle abutting the counter-clockwise edge of the sector (triangle eBC in Fig. |4]) 
with area no more than ^-cot(if). In fact if, for an arbitrary small e, either 5 > sin(3e7r) or 
9 > ip + rjn + 2en, then the area of the intersection between 77-disk and the network region 
is at least (rj/2 + e)nR 2 . Otherwise, it is at least r]irR 2 /2. Thus, after averaging over 5 and 9 
(note that the number of nodes approximately 8R from the edge is uniformly distributed with 
an expected value proportional to 1 — \fd8 when d is very small), the probability that node e is 
disconnected from the rest of the network is derived to be no more than 



(1 - 12ne 2 ) [(1 - 277)e-" diV + 2 V e-^ 2+ ^ dN ] + ^vreV^ 2 < e -^ 2+ ^ dN + 12ne 2 e^ dN / 2 , (2) 

where 127re 2 is an upper bound on the probability that 5 < sin(3e7r) and 9 < + r]ir + 2en, 
and (1 — 2r]) is the probability that the whole empty //-disk falls within the network region. 
Choosing e = 2 l °^ N ^ maximizes the right-hand side of ©, which leads to the upper bound 
48?f l °(dN)^ +1 e~ vdN / 2 for the probability that node e is disconnected from the rest of the network|f|. 
Subsequently, accounting for the expected dN nodes in node e's transmission range and the 
expected 2\fdN edge nodes, we obtain the upper bound on the probability that some edge 
nodes are not connected to the rest of the network as: 

a% < 967Tlog^iV) + l e -^ /2 + 2V - dNe - dN/2j (3) 

where the second summand on the right-hand side is the probability that some edge nodes have 
no other nodes within their transmission ranges. 

C. Calculation of 

Finally, summing (OQ) and ©, we obtain the bound a N on the probability that the network is 
not connected: 

a N < N(dN + l)e^ dN + — ± 1 e -vd*/* + 2 VdNe~ dN / 2 . (4) 

yd 



5 This upper bound is looser but, independent of e. 
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This bound on a N is valid for all N and d = o (iV" 2 / 3 ) , and is asymptotic to 967rl °g( dJV ) e -nWI* 
which goes to zero if dN — > oo as N — > oo. Particularly, assuming d = clo ^ N ^ and r/c > 1 
results in the network being connected with probability approaching one as N — > oo, which 
shows the consistency between our result and the ones derived in [0 and (H for 77 = 1. 

IV. Hop-Count Between S-D Pairs 

Assuming that the network is connected as stated in Definition [Q we now investigate the 
question of how many hops it takes for a packet to traverse from the source to its destination 
with the random 1/2-disk routing scheme. To answer this question we need to have a mathe- 
matical description of the process that models the distance evolution from a typical packet to its 
destination. 

Specifically, we set the intended destination node at the origin and assume that the packet 
starts at the source node X = (—h, 0), where X t is the (Cartesian) coordinate of the packet and 
r t := \\X t \\ is the (Euclidean) distance from the packet to its destination at time t £ IN. 

The packet starts at the source node X with |RSR D that is a 1/2-disk centered at X and 
oriented towards the destination at (0,0). The next relay X\ is selected at random from those 
contained in D (the rRSL rule). This induces a new |RSR D±, also a 1/2-disk but centered at 
X 1 and oriented towards the destination, such that D\ is both a translated and a rotated version 
of D . Relay X 2 is selected randomly among the nodes in D\, and the process continues in the 
same manner until the destination is within transmission range. We claim that the packet reaches 
its destination whenever it enters the transmission/reception range of the final destination, i.e., 
r T < R, for some r £ IN. 

Let \C\ denote the area of a region C, and C c := A — C denote the complement of C with 
respect to network region A. Each |RSR D t has a radius R, with \D t \ = | R 2 . Define St ■— h—r t 
and Y t := St+i — S t , and let N(D t ) be the number of nodes in D t . We would like to know 
how similar Y t (and consequently S t ) is to a Markov process. Note that even though new relays 
are chosen uniformly within each |RSR, the increments Y ,Yx,--- are not independent due to 
the fact that the orientations of all |RSRs are pointing to a common node (destination). This 
results in the dependence of the node distribution in D t , not only on X t , but also possibly on 
Xi, < i < t, where this dependence increases as the packet gets closer to the destination. Note 
that throughout this section we assume that the network is connected, i.e., N(D t ) > for all t, 
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unless otherwise stated. 

Since tracking the dependence of Y t on all JQ, i < t is tedious, in this work we only investigate 
the dependence of Y t on X t and X t _i and in particular, determine how close is the distribution 
of the nodes in D t to a uniform distribution when just knowing the location of the previous 
node X t _i []. Note that conditional on X t (or equivalently on D t ) and network being connected 
(N(D t ) > 0), N(D t ) has a Poisson distribution with intensity A|A|> given it is positive (non- 
zero). What is less clear, however, is the nature of N(D t ), given Xt-i. 

Observe that D t only depends on X t . Given X t , X t ~\, N(D t _ 1 ), and N(D t ) > 0, the number 
of nodes N' t in D t H A-i is Binomial \N(D t -i) — 1, ^^~j^ J + I{x t _ie-D t } an ^ independent 
of the number of nodes N" in D t D A-i' which is Poisson(A|A n -DJLJ); where represents 
the indicator function, i.e., I{e} = 1 if event E happens and I{ E y = otherwise. Moreover, 
conditioned additionally on the two random variables N' t and N", each collection of nodes 
(located in D t fl D t -\ and D t H is uniformly distributed on the respective areas. This 

does not, however, imply that the combined collection of nodes is uniformly distributed on D t . 
Nonetheless, according to the following proposition, the error that results from proceeding as if 
X t+ i is chosen uniformly on D t is negligible for large A. Essentially, knowing X t , the distribution 
of nodes in D t is independent of the location of the previous nodes (Xj, j < t) for large A. 

Proposition 1. Given the location of previous relaying node (X t _i), the distribution of the nodes 
located in the ^RSR of the current node (D t ) approaches the uniform distribution as A — > oo. 
In particular, the conditional probability of selecting the next node from D t fl D t ~i satisfies: 

D t n D t . 



l \ |AnA-i| ( K 
l — ; t L - L , — < E 1 



A|A-i|y |A| " V N(D 



iV(A-i) > 0,N(D t ) > 0,X t ^,X t < 



1 1 



|A| 
(5) 



Proof: Refer to Appendix [A] ■ 
Accordingly, we can model the distance from a packet to its intended destination, {r t } tG]N , as 
a Markov process solely characterized by its step size £ t := —Y t : 

6 The analysis gets more complicated as we consider longer history of the previous relaying nodes, i.e., X t -2, X t -3, 
but we conjecture that the error resulted from neglecting Xt-2, Xt-3, • • • should remain of order O (1/A), where 1/A is the 
error resulting from neglecting Xt-i, as shown in Proposition Q] This is because there are only finitely many previous relays 
whose |RSR might intersect with Dt, and we can assume, as a worst case scenario, that all of them have roughly similar and 
independent effects on the node distribution in D t . 
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r t+1 = r t + £ t I {rt>Rh (6) 

with r = h being the initial distance between the source and the destination. Here we assumed 
that the step size at each time is merely a function of the current position of the packet and the 
effect of the knowledge on the previous locations is negligible, i.e., Pr(£ t < x \ r u ,u < t) ~ 
Pr (6 < x | r t ), which is justified by Proposition \T\ for large A. Fig. |5] depicts (a realization of) 
the evolution of distance of a packet to its destination at time t based on ©. 

The conditional distribution of the step size F(x,r t ) := Pr(£ t < x \ r t ) for —R < x < 
— r t can be obtained as follows. Given r t , the probability that the next hop is at least 
x distance farther/closer to the destination than the current node is equivalent to the probability 
that the relay falls in the intersection of D t and a circle centered at the destination with radius 
r t + x, i.e., the dashed area depicted in Figs. |6fa) and |6{b) for x > and x < 0, respectively. 
Any node in the interior of a circle centered at the destination with radius r t + x is at least 
x farther/closer to the destination than the current transmitting node (which is located at r t ). 
Consequently, we can determine the conditional probability distribution of the packet's next 
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Fig. 6. Derivation of step size distribution (a) for x > (b) for x < 0. 



displacement at time t as: 



F(x,r t ) := Pr (£ t < x|r f ' 



7Ti? 2 



+ m + cc) cos' 



2r t R 

_ x ( 2r t (r t +x) - (R 2 -x 2 ) 



2r t (r t + x) 
-^{R 2 - x 2 ) (4r t (r t + x) - {R 2 - x 2 )) 



- I {x . >0 } ([r t + x) 2 cos 1 ^ - qj— j + r t ^J (r t + x) 2 - r^j 



(7) 



for — i? < a; < i/r 2 + i? 2 — r t . Observe that this distribution directly depends on the distance 
between the packet and the destination in the current time slot, which makes direct analysis of 
the random process extremely tedious. Hence, as shown in Fig. [71 we rewrite © as 



n+i = (n - x' t ) 2 + y' t 2 , 



(8) 



where (x' t , y' t ) is the projection of X t +\ — X t onto the local Cartesian coordinates with the current 
node as the origin and the x axis pointing from the current node to the destination. Based on 
Proposition [Q (x' t , y' t ) is distributed uniformly on D t which makes ([8]) the Markov approximation 
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Fig. 7. Distance between the next relay and the current node projected onto to the local coordinates at the current node. 



Of©. 

Define v r : = inf{£ : r t < r}, for r > R as the first time that the packet gets closer than 
r to the destination. Hence, vr represents the first time the packet enters the reception range 
of the destination and vr + 1 indicates the number of hops that a packet traverses to reach its 
destination. It is easy to show that v r is a stopping time |[T0l . It is clear that 

r — R < r Vr < r. 

Furthermore, let g(r, x', y') = \J (r — x 1 ) 2 + y' 2 — r. We observe that g is a decreasing function 
over r, for fixed (V, y'). Thus, for t + 1 < u r , we have r t > r and 

- x' t < n+i ~r t < g(r, x' t , y' t ). 

Hence, we have 

r-R<r Ur <h + Y,g(r,x't,y' t ), (9a) 
t=i 

I'r 

h + ^2(-x t ) < r Vr < r. (9b) 

t=\ 

Note, as well, that (refer to Appendix |B]) 
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4 R /? 

3^ = E(-x t ) < E(g(r, x' t , y' t )) < E(g(R, x' t , y' t )) < -- < 0. (10) 



Now, applying Wald's equality to ( |9a| ) and (|9b|) and rearranging, we obtain a bound on the 
expected value of the stopping time v r \ 

37r(h-r) .,. h-r + R h Ah 

i '-<E{y r \h)< —, ; K r< —, ; 77T < — • (11) 

4i? ~~ -#(<7(r,a^)) " -£(<7(rX Z/0) " i2 

Substituting r with i? we obtain a general bound for the expected number of hops that a packet 
traverses to reach the destination (minus one) as: 

This implies that the packet reaches its destination almost surely (a.s.) in a finite number of 
steps when the network is connected, which happens with probability no less than 1 — a N as 
obtained in ©. Specifically, with dN — > oo as N — > oo, we have Pr(^ fi < oo) — > 1. 

Moreover, when the ratio h/R (i.e., the ratio between the source-destination distance and the 
transmission range) is large, we can obtain a tighter bound on the expected number of hops be- 
tween a source-destination pair that are h apart. Since r Ur < r, we must have 
E(ur I h) < E{y r | h) + E(ur \ r). Thus, by CCD) and proper substitutions, we have 

gfrg » V < ? + ^ (13) 

for all r G (i?, h). From this, it is easy to deduce that 

hm ; < mm — ; — — = — . (14) 

h/r-^oo h/R r>R —E(g(r,x' t ,y' t )) 4 

This implies that as h/R — > oo, the expected number of hops a packet traverses between a 
source and its destination is asymptotic to 

Furthermore, we can determine the expected initial distance h between an arbitrary source- 
destination pair in the network by deriving the expected distance between any two points that 
are randomly located inside the network area. The problem of quantifying h is well studied in 
the literature [9], with the following known results for two special convex network regions: If 
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the region is a planar disc with diameter a, we have h = 64a/(457r) »s 0.4527a; and if it is a 
square with side length a, we have h — (2 + \/2 + 51og(\/2 + 1)) a/15 « 0.5214a. 

V. Conclusion 

In this paper, we determined the sufficient conditions for a localized routing scheme, namely 
the random 1/2-disk routing scheme, to be convergent based on a convergence notion that 
not only considers the feasibility of packet delivery to the destination via possibly multi-hop 
relaying, but also requires the packet delivery to occur in finitely many hops. To this end, we 
showed that, given certain conditions, there exists a sequence of relays complying with the 1/2- 
disk routing scheme between any two nodes in the network with probability approaching one 
as the (network) node density goes to infinity. Furthermore, we showed that given a network 
satisfying those conditions, the random 1/2-disk routing scheme will deliver the packets to their 
destinations in finitely many hops almost surely, by modeling the distance evolution from a 
packet to its destination as a Markov process. 

Appendix A 
Proof of Proposition CD 

First, let us consider the distribution of a Poisson point process conditioned on deleting one 
point. Let N be a homogeneous Poisson point process with intensity A and fix a region D. 
If N(D) > 0, one point in D is selected at random and removed. Let X be the location of 
that point. The distribution of N on D c remains Poisson and independent of N on D, and thus 
independent of both N(D) and X. Let N' be the (point) process with the point at X deleted. 

Let Ai,A 2 , ■ ■ ■ , A k be a partition of D. Given N(D) > 0, the points in D are distributed 
uniformly. If one point is removed at random, the remaining points are still distributed uniformly 
on D. Hence, 



Pr\^]{N'(A J ) = n J } | N(D)>0,X 



i=l j=l 

X\D\ 



(\\A 3 \) 



nj+l {j=i} 



Ai\ n i + Mi =i > 



-A|^| 



77 



+ I{i=i})! 



ri\ + • • ■ + n k + 1 



l_ e -A|D|)( ni + ... +nfc 



nil 



n 



'3 J- 



(15) 
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since \Ai\ + • — \-\Ak\ = \D\. Therefore, conditional on N(D) > 0, N' is independent of X. In 
particular, 



Pr (N'(D) 



n 



N{D) > Q,X 



(MD\r +1 



(n + l)!(l 
Pr (N(D) 



-A|£>p 



n 



-X\D\ 



N{D) > 



Furthermore, given n-y + ■ ■ • + n k = n > 0, we have 



Pr ( f]{N'(Aj) = rij} N(D) > 0, N'(D) =n,X 



\j=i 



n 

nx---n k 



n 

i=i 



Thus, for A C D and given iV'(D) = n > 0, iV(A) is conditionally binomial (n, Without 
knowing N'(D) > 0, however, we obtain from <TT3T > that 



Pr (a? 7 (A) = A; iV(D)>0,X 



E 



1 3=0 J 



A|D|e- A l D l (\\A\Y 



fc! 



(16) 



where the second equality is due to 




After the aforementioned preliminaries, we now proceed with the proof of Proposition [T] 
Suppose C is a random set that depends only on X. The points of N', if any, that are in 
CD := C fl D are uniformly distributed and independent of the points in CD C , which also are 
uniformly distributed (if any). The combined points are uniformly distributed on C only if the 
expected proportion of points in CD is 

However, it is not in our case as we now compute. Given N'(C) > 0, the probability that 
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a randomly selected point in C also is in D is E(^^ | N'(C) > 0,N(D) > 0,X). Let 



N'(CD) 
N'(C) 



when N'(C) = 0. We have 



Pr(N'(C)>0 N(D) > 0,x) = 1 -Vt(n'(CD) = 0,N'(CD c ) = N(D)>0,X > ) 

_ X\D\e-^e^\-l _ A|OTC| 
l-e- A l D l A|C C L>| 



|D| 1-e 



-A|C C D| 



|C C D| 1 — e -^l^l 



-A|C| 



and so 



\C C D\ 



" A|c| < Pr (N'(C) > N(D) > 0, < 1 



-A|C| 



(17) 



Using the observation above and (1161) we have, 



N'(C D) 

E( „ V ,„/ V(c)>o> | N(D) > 0,X) 



N'(C) 



Therefore, 



y n e X\C^D\y dy iM£^r e -\\CD^\ 



\\D\e- x W ~ ~ n (\\CD\) n 
1 - e~ A l D l n _|_ m n i 

n=0 m=0 

n=l ' 

A | jD | e -A|CUD| 

1 - e" A l D l 
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\CD\y 



\\D\ e -*\CUD\ rl 

1-e-m Jo \CD\y+\CDt 



^ e \(\CD\y+\CD-\) _ -A e A|C c D|y d2/ 



| CD | A^Dje A I CUD I r 1 , X (\CD\y+\CD°\) _ ,\ \\C-D\yi 

~ \C\ l-e^l J [e 1)6 ay 

\CD\ 



\C\ 



Pr(V(C)>0 N(D)>0,X). 



(18) 



ESl^)>0,^)>0,X)< iCT| 



iV'(C) 



|C| 



Noting that 
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ae~ a / (l-y)e a ydy= <-, 

Jo a a 

we also may ascertain that 

\CD\ 



\c\ 



So, if A is very large, the selected point is slightly less likely to be in D than would the case if 
we could assume N' is Poisson on C (or that the points in C truly were uniformly distributed). 

Appendix B 
Derivation of inequality (flOl ) 

We have (x' t , y' t ) = (Rv cos(6 ) ), Rv sin(6 1 )), where 9 ~ uniform(— ir/2, 7r/2) and v ~ beta(2, 1) 
are independent. Thus, 

2 r^ 2 r 1 ar 

E(x' t ) = R- / cos(#) d9 / 2v 2 dv = — . 

Jo Jo 3tt 

Also, by first changing x to 1 — x and then using polar coordinates, 

E(n(R.x',.v'.)) + 1 = - / 

R 



^E(g(R,x' t ,y' t )) + l = - [ [ \ x ^< x y/(l - xf + y*dxdy 

ft 77 Jo Jo 

= - [ I l (1 _ a .)2 +y 2< 1 \Jx 2 + y 2 dxdy 

71 Jo Jo 



- / / 2t> 2 dwdfl + - / / 2t? 2 dwdfl 

K Jo Jo 71 Jtt/4 Jo 

3^ 



y (jsecc?) 3 + (2sin#) 3 )dfl 
3(2 3 / 2 ) + 61og(l + y/{2)) + 64 - 5(2 7 / 2 ) 



0.7499728. 



9tt 

Hence ^ y' t )) < -f . 
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